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Theoretical determination of 775's electromagnetic decay width 

Nicola Fabiano a 

a Pcrugia University and INFN, via A. Pascoli 1-06100 Perugia, Italy 

We discuss the theoretical predictions for the two photon decay width of the pseudoscalar etab meson. Pre- 
dictions from potential models are examined. It is found that various models are in good agreement with each 
other. Results for etab are also compared with those from Upsilon data through the NRQCD procedure. 
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1. Introduction 

First evidence of r/t, was given in 2001 by 
ALEPH . The search for ijb is done in the decay 
channel into two photons. One candidate event is 
found in the six-charged particle final state and 
none in the four-charged particle final state, giv- 
ing the upper limits: T(rjb — > 77) x BR(rjb — > 
4 charged particles) < 48 eV , r(?7b — > 77) x 
BR(jjb — + 6 charged particles) < 132 eV . This 
observation was the main motivation for a the- 
oretical estimate of electromagnetic decay width 
for r)b [23- More recently CDF has done a search 
for 77b in the process [3] rjb — > J/tjiJ/tj) with 7 
events where 1.8 events are expected from back- 
ground. 

This note is organised as follows: in Sect. 2 we 
shall compare the two photon decay width with 
the leptonic width of the T. Sect. 3 is devoted 
to the potential model predictions for rjb — > 77, 
with potential given by |5I6I7| . In Sect. 4 we 
show the predictions for 775 decay widths, using 
the procedure introduced in [H] for the descrip- 
tion of mesons made out of two non relativistic 
heavy quarks, by means of the Non Relativistic 
Quantum Chromodynamics NRQCD. In Sect. 5 
we compare these different determinations of the 
775 — > 77 decay width together with a result based 
on a recent two-loop theoretical analysis of the 
charmonium decay 



2. Relation to T electromagnetic decay 
width 

We start with the two photon decay width of a 
pseudoscalar quark-antiquark bound state in the 



singlet picture, generically written as 
T = NP x P 



(1) 



NP being the nonperturbative part and P the 
perturbative correction. The nonperturbative 
part for a state with a given £ takes the form 
NP ~ \d e ip(0)/dr e \ 2 . The wavefunction is ob- 
tained by a solution of a Schrodinger equation 
with a suitable potential V(r). The perturba- 
tive expression is given by P — F(a s (m q )) and 
comes from the on-shell matrix element. The 
two photon decay width of a pseudoscalar quark- 
antiquark bound state |10) is given by 



T B (Vb -►77) = 12QV47T 



M 2 



1 B 



(2) 



M = 2mb + Eb being the mass of the state, Eb the 
(negative) binding energy; with first order QCD 
corrections |11| . which can be written as 



rfa -» 77) = 



1 + 2; 
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A first theoretical estimate for this decay width 
can be obtained by comparing eq. ® with the 
expression for the electromagnetic decay for the 
vector state T i.e. 



T B (T 



= 4Q 2 a 2 47r 
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and the one-loop complete formula 
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Using the expressions in eqs. © and JSJ we can 
estimate the rjb decay width from the the mea- 
sured values of the leptonic decay width of T. 
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We assume the wavefunctions of the two states 
to be equal, that leads to an error proportional 
to (a s /m\) PI- Taking the ratio of the eqs. © 
and © and expanding to first order in a s , we 
obtain: 



r(Tk^yy) 



T{ Vb 77) 1 (1 - 3.38a s /7r) 



r(T^e+e-) 3 (1 - 5.34a s/ V) 



1 



1.96— + 0(c? s 
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(6) 



In order to compute the correction we start from 
the two-loop expression for a s and the value at 
the Z mass [T5] a s {M z ) = 0.118 ± 0.003. Using 
the renormalization group equation to evaluate 
a s {Q = 2m b = 10.0 GeV) = 0.178 ± 0.007, and 
the latest measurement of the T 



r exp (T 



e e 



■) = 1.32 ±0.05 keV 



(7) 



for the r\h one obtains 



r(r? 6 -> 77)±Ar(7y 6 -> 77) = 489±19±2 eV ,(8) 

where the first error comes from the uncertainty 
on the T experimental width, the second error 
from a s . We have assumed the a s scale to be 
Q = 2mb = 10.0 GeV. This choice is by no way 
unique as shown for the r\ c decay and in fig. Q 
we show the dependence of the r\h photonic width, 
evaluated from eq. |7|. upon different values of 
the scale chosen for a s . 

Unlike the r\ c case there are no experimen- 
tal measurements of this decay; we shall assume 
therefore that it is not possible to determine a 
scale choice of a s for the % decay. We will have 
to include this fluctuation in the indetermination 
due to radiative corrections. 

3. Potential Model predictions 

We present now results for rjb — > 77 from the 
potential models. This method allows us to ob- 
tain the absolute width, through the wave func- 
tion at the origin computed from a Schrodinger 
equation. For the calculation of the wavefunc- 
tion JHj we have used four different potential 
models, starting from the one of Rosner et al. |S] 
V(r) = X((r/r a ) a - l)/a + C . with r = 




Figure 1. The dependence of the r\b decay width 
to 77 (in eV) is shown with respect to the scale 
chosen for a s in the radiative corrections. The 
fluctuation is of order 2% 



— 1.305 GeV, and the QCD inspired potential Vj 
of Igi-Ono |6TT7| 

Vj(r) = V AR (r) + dre- gr + or, 
4ai 2) (r) 



V AR {r) 



(9) 



1 GeV" 



= -0.14 , A = 0.808 GeV, C = 



with two different parameter sets, corresponding 
to = 0.5 GeV and = 0.2 GeV re- 

spectively 6 . We show also the results from 
a Coulombic type potential with the QCD cou- 
pling a s frozen to a value of r which corresponds 
to the Bohr radius of the quarkonium system, 
re = 3/(2mi,a s (rB)) (see for instance jjj). The 
latter has the advantage of providing analytical 
solutions for the wavefunction and the energy 
levels: E n = -4ma s (r B ) 2 /(9n 2 ) and |V>(0)| 2 = 
(2ma s (r B )/3) 3 /7r 

We have to stress the fact that the scale of a s 
occurring in the radiative correction and the one 
of Coulombic potential are different. We show in 
fig. (J5J) the predictions for the decay width from 
these potential models with the correction from 
eq. © at an a s scale Q = 2m,b . 

For any given model, sources of error in this 
calculation arise from the choice of scale in the ra- 
diative correction factor discussed before and the 
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■ Rosner et al. 
Igi-Ono, A MS = 0.2 

. Igi-Ono, A MS =0.5 

■ Coulomb 



that is a sum of terms in inverse powers of m g , 
each of which factorises into a perturbative coef- 
ficient 3/„(a s ) and nonperturbative matrix ele- 
ment (O n )v,p- 

According to BBL, in the octet model 
for quarkonium, the electromagnetic and light 
hadrons (LH) decay widths of bottomonium 
states are given by: 



5.0 

Bottom mass (GeV) 



T(T -> LH) 



2(T|Oi( 3 5i)|T) fl0_ 2 _ HT 
243 ^ 27, 



Figure 2. The dependence of r\ b decay width to 
77 in eV for different potential models is shown 
as a function of mj. The perturbative correction 
scale is here set to be Q — 2m b . 
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2<T|P 1 ( 3 S 1 )|T) 17.32 x [20(^ 2 - 9)] 3 



choice of the parameters. The difference of the 
prediction between the various models introduce 
another indetermination for the decay value. We 
can estimate a range of values for the potential 
model predictions for the radiative decay width 
r (% -> 77), namely: 

r(ry fc -> 77) = 466 ± 101 eV , (10) 
that is an indetermination of around 22%. 

4. Octet component procedure 

We will present now another approach which 
admits other components to the meson decay be- 
yond the one from the colour singlet picture (Bod- 
win, Braaten and Lepage) 0. NRQCD has been 
used to separate the short distance scale of anni- 
hilation from the nonperturbative contributions 
of long distance scale. This model has been suc- 
cessfully used to explain the larger than expected 
J/4 1 production at the Tevatron and LEP. The 
decay width expression for a pseudoscalar P or a 
vector state V is given by means of NRQCD from 
the expansion: 



V,P 



E 



25sf n (a s 

dn-4 



■(O r , 



>V.P 



(11) 



r(T - 

13 as 

17 



e+e") 
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ira 2 (14) 



2{ Vb \O l ( 1 S )\m) 



r(% -> 77) = " vl ° rn 2 "" m ° / fQV [i+ 

2(r lb \P 1 ( 1 S )\r lb )A . 2 
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IT 
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There are four unknown long distance coeffi- 
cients, which can be reduced to two by means of 
the vacuum saturation approximation: 

Gi = (TldfSOIT) = (rjblOxCSo^rib) (16) 
F x = (TlPi^SiJIT) - (vblPxCSo^Vb) (17) 

which is correct up to 0(v 2 ), where v is the quark 
velocity inside the meson. In the potential model 



4 



language this translates to consider the two wave- 
functions of the pseudoscalar and the vector state 
to be equal. With this position we obtain a sys- 
tem of equations for the nonperturbative coeffi- 
cients Gi and F\ of T electromagnetic decay and 
light hadrons widths which in turn allows us to 
compute the rjb decay widths. 

The BBL approach gives the following decay 
widths of the r\b meson: 

r(??fc -> 77) = 364 ± 8 ± 13 eV (18) 
and 

T(rib -> LH) = 57.9 ± 4.6 ± 2.8 keV , (19) 

where the first error comes from the uncertainty 
on the T experimental width, the second error 
from a s . The improvement of the error on 
eq. ifH^I with respect to the previous analogous 
determination on the r\ c decay @] is due to better 
error on the experimental measures of the T de- 
cay widths compared to the one of the J/ip, and 
the smaller indetermination on the a s value due 
to the higher energy scale involved in the decay. 
These reasons, together with the fact that the 
potential models used are fitted for the cc sys- 
tem, justifies the improvement of accuracy given 
in eq. I|18|l compared to the one of eq. I|1U|) . 

5. Comparison between models 

For comparison we present in fig. @ a set of 
predictions coming from different methods. 

Starting with potential models, we see that the 
results are in good agreement with each other. 
The advantage of this method is that we are giv- 
ing a prediction from first principles, without us- 
ing any experimental input. Since there are cur- 
rently no experimental measures for the rjb — > 77 
decay, we shall use this prediction as a reference 
point, as it has proven to be reliable in the case 
of charmonium decay @]. The second evaluation, 
given by BBL using the experimental values of 
the T decay, is on the left limit of the potential 
models value. This is true also for the determi- 
nation of the BBL procedure with nonperturba- 
tive long distance terms taken from from the lat- 
tice calculation [T^J, affected from a large error of 
around 30%. The advantage of the latter is that 
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Figure 3. The rjb two photon width as calcu- 
lated in this paper using (starting from below) 
Potential Models results, BBL procedure with in- 
put from J ftp decay data, Lattice evaluation of 
G\ and F 1 factors, Singlet picture with G\ ob- 
tained from T — > e + e~ and T — * LH processes 
respectively, and the two-loop enhanced proce- 
dure. The legend shows the percentage error on 
each value. 



its prediction, like the one from potential mod- 
els, does not make use of any experimental value. 
Next is the point given by the singlet picture from 
the electromagnetic decay of the T, aligned with 
the aforementioned results of the BBL procedure. 
The point above is obtained also from the singlet 
picture with the T decay into light hadrons, in 
agreement with the results given from the poten- 
tial models. We notice that in analogy to the 
charmonium case (see 0] and references therein) 
the singlet results obtained from the T decay are 
in disagreement with each other, in this case by 
only la. The last point from a two-loop enhanced 
calculation given by |9I1| is in agreement with the 
potential model result and the singlet decay from 
the T — > LH process. 

6. Conclusions 

The r(?y/ J — > 77) decay width prediction of 
the potential models considered gives the value 
466 ± 101 eV, in agreement with the naive esti- 



5 



mate from the T decay given by ©• Predictions 
of the BBL procedure are consistent with the po- 
tential model results, for both the long distance 
terms G\ and F\ extracted from the T experimen- 
tal decay widths and the one evaluated from lat- 
tice calculations. The results from the singlet pic- 
ture are also consistent with the potential model 
results. Finally the two-loop enhanced prediction 
is in good agreement with the potential model re- 
sults. 
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